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Current-driven instabilities in weakly ionized disks 
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Cool weakly ionized gaseous rotating disk form the basis for many models in astrophysics objects. Instabilities against 
perturbations in such disks play an important role in the theory of the formation of stars and planets. Traditionally, axisym- 
metric magnetohydrodynamic (MHD) and recently Hall-MHD instabilities have been thoroughly studied as providers of 
an efficient mechanism for radial transfer of angular momentum, and of density radial stratification. In the current work, 
the Hall instability against axisymmetric perturbations in incompressible rotating fluid in external poloidal and toroidal 
magnetic field is considered. 
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1 Introduction 

The origin of turbulence in astrophysical disks is often at- 
tributed to hydrodynamic and hydromagnetic instabilities 
that can occur in differentially rotating stratified gas. The 
magnetorotational instability (MRI) is usually considered 
as one of the possible candidates because it can operate in 
a conductive flow if the angular velocity decreases with the 
cylindrical radius (Velikhov 1959; Chandrasekhar 1960). 
The MRI has been studied in detail for both stellar and 
accretion disk conditions (see, e.g., Fricke 1969;Safronov 
1969;Acheson 1978,1979;Balbus & Hawley 1991; Kaisig, 
Tajima & Lovelace 1992;Zhang, Diamond & Vishniac 1994 
Note that the MRI occurs only in the presence of a weak 
magnetic field, because a sufficiently strong field can sup- 
press the instability completely. Numerical simulations of 
the MRI in accretion disks ( Hawley, Gammie & Balbus 
1995; Brandenburg et al. 1995; Matsumoto & Tajima 1995; 
Torkelsson et al. 1996;Arlt & Rudiger 2001) show that the 
turbulence generated can enhance essentially the angular 
momentum transport. 

The properties of turbulence in low conductive proto- 
stellar disks can differ essentially from those of accretion 
disks. The magnetic Reynolds number is not very large in 
protostellar disks and, hence, the field cannot be treated as 
'frozen' into the gas ( Gammie 1996). The influence of 
ohmic dissipation on the MRI has been considered in the 
linear (Jin 1996) and nonlinear regimes ( Sano Inutsuka & 
Miyama 1998). Calculations of many authors indicate that 
the MRI is unlikely to be the source of turbulence in pro- 
tostellar disks since it arises only in a highly conductive 
plasma. For instance, Turner, Sano & Dziourkevitch (2006) 
have considered in detail turbulent mixing caused by the 
MRI in protostellar disks and argued that the MRI does not 



arise in the midplane even under the most favorable con- 
ditions because the midplane is shielded from cosmic rays 
which are the main ionizing factor. However, the number 
of instabilities that occur in disks is not restricted by the 
MRI alone. An analysis of MHD modes in stratified disks 
demonstrates a wide variety of instabilities even in the case 
of a very simple magnetic geometry ( Keppens, Casse & 
Goedbloed 2002). More complex magnetic geometries with 
a non-vanishing radial field lead to additional instabilities 
( Bonanno & Urpin 2007, 2008). Generally, even a pure hy- 
drodynamic origin of turbulence cannot be excluded (see, 
e.g., Dubrulle et al. 2005; Lesur & Longaretti 2005). 

It was first pointed out by Wardle (1999) that poorly 
conducting protostellar disks can be strongly magnetized if 
electrons are the main charge carriers. As a result, trans- 
port must be anisotropic with substantially different prop- 
erties along and across the magnetic field. The effect of the 
magnetic field on transport properties of plasma is charac- 
terized by the magnetization parameter a c = ojbt where 
ujb = eB/m c c is the gyrofrequency of electrons and r is 
their relaxation time (see, e.g., Spitzer 1978). In protostellar 
disks, t is determined by the scattering of electrons on neu- 
trals and can be calculated by making use of the fitting ex- 
pression for the cross-section obtained by Draine Roberge 
& Dalgarno (1983). Then, we have for the magnetization 
parameter 

-1/2 



a c « 21 Bn 14 T 2 



(1) 
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where B is the magnetic field measured in Gauss, ni4 = 
n/10 14 cm" 3 and T 2 = T/100 K with n and T being the 
number density of neutrals and temperature, respectively. 
If a c > 1, i.e. B > 0.048 n U y/% G, then the electron 
transport is anisotropic, and the magnetic diffusivity is rep- 
resented by a tensor. In a weakly ionized plasma of proto- 
stellar disks, the difference between the components of the 
magnetic diffusivity which are parallel and perpendicular 
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to the magnetic field is small (see, e.g., Balbus & Terquem 
2001), but the Hall component that is perpendicular to the 
both magnetic field and electric current can be much greater. 
The Hall component of diffusivity is given by a c rj where 
j] = c 2 m e /4:ire 2 n e T is the magnetic diffusivity at B = 0; 
m e and n e are the mass and number density of electrons, 
respectively. Using the fit for the cross-section by Draine 
Roberge & Dalgarno (1983), we obtain for the magnetic dif- 
fusivity 

T) = 2.34 x 10 3 x- 1 T^ /2 cm 2 s _1 , (2) 

where x e = n e /n is the ionization fraction. 

A stability analysis of the MRI done by Wardle (1999) 
shows that the Hall effect can provide either stabilizing or 
destabilizing influence depending on the direction of the 
field. A more general consideration of the Hall MRI has 
been done by Balbus & Terquem (2001). They found that 
the Hall effect changes qualitatively the stability properties 
of protostellar disks and can lead to instability even if the 
angular velocity increases outward. These authors, however, 
did not take into account the effect of gravity that is cru- 
cial for disks. A consistent consideration of the linear MRI 
under the combined influence of the Hall effect and grav- 
ity has been done by Urpin & Rudiger (2005) who derived 
also the criteria of several other instabilities that can occur 
in protostellar disks. The properties of the MRI modified 
by the Hall effect has been considered also by Salmeron & 
Wardle (2005). These authors argued that the MRI is active 
in protoplanetary disks over a wide range of field strengths 
and fluid conditions. The Hall conductivity results in a faster 
growth of perturbations and extends the region of instabil- 
ity. Recently, Liverts Mond & Chernin (2007) and Shtem- 
ler, Mond & Liverts (2007) have considered the Hall in- 
stability (HI) for non-axisymmetric perturbations. This in- 
stability differs from the MRI and it results from the fast 
magnetosonic waves in contrast to the Alfven nature of the 
MRI. The HI instability is proposed as a viable mechanism 
for the azimuthal fragmentation of the protoplanetary disks 
and planet formation. The non-axisymmetric instability is 
caused basically by the combined effect of the radial strati- 
fication and Hall electric field. 

Apart from the Hall effect, the stability properties of 
magnetic protostellar disks can be influenced by a num- 
ber of other factors, for example, the electric currents. In 
differentially rotating disks, the azimuthal field is gener- 
ated by stretching from the poloidal one, and poloidal cur- 
rents are necessary to maintain this azimuthal field that typ- 
ically is not current-free. The generated azimuthal field can 
be stronger than the poloidal one if the magnetic Reynolds 
number is larger than 1 . The azimuthal field and associated 
currents can be important for stability of disks even if the 
Hall effect is negligible (see Pessah & Psaltis 2005). In pro- 
tostellar disks, the effect of currents maintaining the mag- 
netic configuration is accompanied often by the Hall effect 
that changes crucially the stability properties. 

In the present paper we consider the linear stability 
properties of magnetic protostellar disks taking into account 



the combined influence of the Hall effect and electric cur- 
rents. The criteria of instability are derived, and the growth 
rate of the various modes is calculated. 

2 Basic equations 

Consider the stability of a magnetized protostellar disk of 
a finite vertical extent. For the sake of simplicity, the un- 
perturbed angular velocity is assumed to be dependent on 
the cylindrical radius s alone such as £1 = f2(s); (s, <p, 
z) are cylindrical coordinates. The magnetic field, B = 
(B s , B v , B z ), is assumed to be weak in the sense that the 
Alfven speed, ca, is small compared to the sound speed, c s . 
This enables us to employ the Boussinesq approximation for 
a consideration of slowly varying modes. In the unperturbed 
state, the disk is assumed to be in hydrostatic equilibrium in 
the s- and z-directions, 
Vp „ 1 



G 



-rotB x B 



G 



- - 9 • (3) 

p Anp 
here g is the gravity force per unit mass. It should be noticed 
however that the pressure gradient in the unperturbed disks 
is mainly determined by gravity and centrifugal forces. The 
unperturbed Lorentz force namely the second term on the 
r.h.s. Eq. (0 is usually much smaller than gravity and cen- 
trifugal forces (the first term on the r.h.s. of the last equa- 
tion) and therefore may be neglected. This is due to the as- 
sumption of small thickness of the disk. Detailed asymptotic 
description of that fact may be found in Regev (1983);Kluz- 
niak & Kita (2000);Shtemler, Mond & Liverts (2007). 

We consider the stability of axisymmetric short wave- 
length perturbations with space-time dependence exp(7i — 
ik ■ x) where k — (k s ,0,k z ) is the wave vector. The 
linearized momentum and continuity equations read in the 
Boussinesq approximation 

ikp 



jv + 2tt x v + e v s(v ■ V)0 

P 

-L[(B-b)k-(k-B)b]-^^ 



- o.GT x 



1 



-J x b, (4) 

cp 

k-v = 0, (5) 

where v, b, p and T\ are the perturbations of the hydro- 
dynamic velocity, magnetic field, pressure and temperature, 
respectively; a = —(d\np/dT)p is the thermal expansion 
coefficient and e v is the unit vector in the azimuthal di- 
rection. It is assumed in Eq. (0]) that the density perturba- 
tion in the buoyancy force is determined by the temperature 
perturbation alone in accordance with the main idea of the 
Boussinesq approximation, p\ = —paT\, and the unper- 
turbed Lorentz force in Eq. ©is neglected. We took into 
account the effect of electric currents (the last term on the 
r.h.s. of Eq. (01). In short wavelength approximation, this 
term seems to be smaller than the previous one by a fac- 
tor ~ kL 1, where L (~ s) is the length-scale of the 
unperturbed magnetic field. In differentially rotating discs, 
however, the electric current J = (c/47r)V x B is mainly 
determined by the (/^-component of the magnetic field that 
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can be substantially stronger than B s and B z because of 
stretching the magnetic field lines in the azimuthal direction 
by differential rotation. If B v satisfies the condition 

B v > kL max( J B s ,B Z ), (6) 

then the effect of electric currents cannot be neglected in 
Eq. ©. For the sake of simplicity, we assume that the 
toroidal field B v depends on s alone, then 

J = J z e Zl J z = -£-^-( a B v ). (7) 
ins os 

Note that, calculating the perturbation of the electric current 
j = (c/47r)V x b in Eq. (|4), we can neglect terms of the 
order of 1/s and assume j w — (ic/Air)k x b. Using Eq. ©, 
we can calculate p from Eq. Then, we obtain for the 
momentum equation 

2k. 



jv + 2f2 xv — —rk ■(flxn) + e^sil'vg 



-aTi 



k 2 



G) 



B)b 



1 

cp 



Jxb-^k-(Jxb) 

k 



(8) 



It is clearly seen from the last equation that the third term 
on the r.h.s. which results from the unperturbed electric cur- 
rent can be larger under condition © than the second term 
which is usually taken into account in a stability analysis of 
disks and which contains only poloidal components B s and 
B z of the background field. 

Since the thermal conductivity of protostellar disks is 
low because of a low temperature (T ~ 10 — 10 3 K), we 
adopt the adiabatic equation to describe the evolution of 
temperature perturbations, 

7T1 + v ■ (AWT) = , (9) 

where (AVT) = VT - V ad T is the difference between 
the actual and adiabatic temperature gradients. Substituting 
Ti into Eq. (0, we obtain the equation that contains only 
perturbations of v and b. The s- and (^-components of this 
equation read 



(7 H -)v 8 — 2uVtv^ 



7 

7*V 
where 



Airp 



(k ■ B)b s 



(20 + sfi') 



-qAVT ■ 



47T/9 



(k ■ B)b^ 



cp 

. J, 



cp 



(11) 



k 2 



(fc.G) 



and /i = k 2 /k 2 . 

As it was mentioned, the effect of the magnetic field 
on kinetic properties of plasma is usually characterized by 
the magnetization parameter a c (see Eq. (Q~|i). We consider 
the most interesting case for protoplanetary disks when this 
parameter is moderate, kL ^> a e . Under this assumption, 
the linearized induction equation reads 



where rj is magnetic diffusivity and n e is the number 
density of electrons. For short wavelength perturbations, 
Eq. (TT2l and the divergence-free condition read 



(7 + iuj A )b 



-e z b. 

en. 



iv(k ■ B) — e 
c(k ■ B) 



Alien,, 



d_ 

os 

k x b + e 



B ■ 



ick 



B' w b v , (13) 
k b = 0, (14) 



v ^en e * v 



where lj a = k z J z /en e — ink 2 , B' = dB v /ds, and J z = 
dJ z /ds; the last three terms on the r.h.s. of Eq. ( fT3l ) together 
with the first term in w A represent the Hall effect that can 
be important in protoplanetary disks. Under condition (|6), 
induction equation ( fT3l can be simplified because the fourth 
term on the r.h.s is small compared to the last term and do 
not influence the behavior of perturbations. Then, we have 
for the induction equation 



d (B 



sfl'b s — iujftb, 



(7 + iw A )b = —iv(k ■ B 
c(k ■ B) 



Alien,, 



kxb, (15) 



where ujh = ck z B' v /4iren e . The s- and (^-components of 
this equation are 



■ uj w b v , 



(7 + iuj )b v 



where 



(7 + iu A )b s = —i(k ■ B)v s 

ik-B k 2 

; Sil V s - -7TUJ w b s 

7 + i<jJ A k% 



ck z Bu 



^0 — — 



ink , uj u 



ck z (k ■ B) 



(16) 



(17) 



(18) 



2iren e s Aiten e 

Eqs. ( [Tol l. (fTTT i. ( TTol l. and (fTTT i describe the eigenmodes 
that exist in a magnetized fluid in the presence of electric 
currents. 



fyp, (10) 3 Dispersion equation and stability criteria 



The general dispersion equation for the set of Eqs. ( [Tol l. 
( fTTT i. (fTol l. and (TTTb is rather cumbersome. Therefore, we 
consider only a particular case of perturbations with the 
wavevector perpendicular to the unperturbed magnetic field, 
k ■ B = 0. The MRI does not occur for such perturbations, 
and they should be stable if the Hall effect is neglected and 
J = 0. Therefore, the particular case k ■ B — allows 
to study a destabilizing influence of both the unperturbed 
electric current and the Hall effect in a situation when other 
factors can cause only a stabilizing influence. If k ■ B = 0, 
then the dispersion equation takes the form 



7 



a 2 7 



- ai7 + «o = 0, 



(19) 



7b = -77V x (V x b) + V x (v x B) 

c 



V x (sOe v x b) w here 



47ren. 



■V x [(V x 6) x B + (V x B) x 6], (12) 



a 2 = iuo, 01 = Ldg+p(K 2 +u)j), a = iuj (u>\ 2 + p,K 2 ) .(20) 
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The characteristic frequencies in these expressions are 
given by 



UJq = LL>i — ILU V , U)\ 



ck z B,r, 2 

, uj v = rjk , 

Ziren e s 



k 2 = 2il(2n + sQ'), u 2 j = - — ^-(^A (21) 

cp os \ s J 

For k- B — 0, only three non-trivial modes exist in the flow 
but other modes are degenerate. 



3.1 Stability in the case s e «l 

If the Hall parameter is small, the Hall effect does not influ- 
ence the stability properties. In this case, ujq w —iui n . Then, 
the dispersion equation takes the form 

7 3 + b 2l 2 + &i7 + b o = 0, (22) 
where all coefficients of this equation are real, 

b 2 =UJ r] , bi = fl! = UJ 2 + (l{K 2 + UJ 2 ), 



b a = uj v (uj 2 + (IK 2 ) 



(23) 



The condition that at least one of the roots of Eq. ( 1221 has a 
positive real part (that corresponds to instability) is equiva- 
lent to one of the following inequalities 

b 2 < 0, M 2 < b , b <0 (24) 

being fulfilled (see, e.g., Aleksandrov, Kolmogorov & Lau- 
rentiev 1985). Since u n > 0, the first condition b 2 < will 
never apply. The other two conditions yield 



(loj^uij < , u 71 (uj 2 + (IK 2 ) < 



(25) 



Both criteria are proportional to the dissipative frequency 
(that is positively defined quantity) and appear only if one 
takes into account magnetic diffusivity. Therefore, both cri- 
teria describe diffusive instabilities that can be relatively fast 
in protostellar disks. In ideal magnetohydrodynamics, we 
have uj v = and bo = b 2 = 0, and dispersion relation (fT9l 
transforms into 

7 2 + ax = 0. (26) 
The condition of instability reads a\ < 0, or 

u 2 g + »(K 2 +uj)<0, (27) 

that differ substantially from dissipative criteria (1251 1. For in- 
stance, to satisfy condition d27l i in disks, energy of the mag- 
netic field should be comparable to the rotational or gravita- 
tional energy of gas. The first condition (EB! can be fulfilled 
even in a much weaker field. Since u n > and (i > 0, we 
obtain from Eq. d25l l the following conditions of instability 



ujj<0, 



UJ 2 + (IK 2 < 0. 



(28) 

■ , (29) 

Eq. d2"9l is the standard criterion of convection modified by 
rotation and usually is not satisfied in astrophysical discs. 
Eq. (l28l is the condition of an instability that can occur due 
to the presence of electric currents (see also Velikhov 1959). 
This instability is associated only with the distribution of 
electric currents in disks. Condition d28l can be rewritten as 



B% - s 2 B' 2 < 0. 



Therefore, the current-driven instability arises if B v de- 
creases with s faster than 1/s or increases outward faster 
than s. Note also that the current-driven instability does not 
depend on the angular velocity profile and can occur for per- 
turbations that are not subject to the MRI since k _L B. 

Since the coefficients of equation (l22l are real there ex- 
ist three real roots or one real and two complex conjugate 
roots. The number of roots with a positive real part is de- 
termined by Routh criterium (DiStefano III, Stubberud & 
Williams 1994), which states that the number of unstable 
modes of a cubic equation d22l is given by the number of 
changes of sign in the sequence 



1, bo 



b 2 bi - bo 



, bo 



(31) 



(32) 



For coefficients (23), this sequence reads 

{1, U! v , UJj, W n {w 2 + (IK 2 )}. 

If the disc is convectively stable (ui 2 + (ik 2 > 0), we obtain 



that under the condition (1281 1 (or ( I30t ) two complex conju- 
gate modes are unstable. If the disc is convectively unstable 
(u 2 + (ik 2 < 0) and the condition (l28l holds then there 
should be only one unstable mode. Only one mode is unsta- 
ble also in the case when the disc is convectively unstable 
but the condition ( f28l ) is not fulfilled. 

The roots 7, (i — 1, 2, 3) of the cubic equation ( [T9| > can 
be represented as ji — xi — a 2 /3. The expressions for Xi 
are 

1 %/3 

Xi=U + V , 2)2,3 = - + V ) ± i ~^ U ~ W ) ' ( 33 ) 

where 



2 

27" 



2? = 7^4 



{u,v) = {- q ±^F^yi\ 

1 „ 1 o 

3p 



a 



ai a, 

1 3 2 



(34) 



(see, e.g., Bronstein & Semendyaev 1957). In a particular 
case a c < 1 when the coefficients of a cubic equation are 
given by Eq. (1231 1. we have 

1 , 1 

V 



(IK 



-(lUJj + -ui v 



1 

P= 3 



[IK 



2 ^-2 



In the limit of small u) n , we have for the roots 



7i 



72,3 = ±i 



to v (ul + (IK 2 ) 



(IK^ 



(iwy 



(UK* 



2uj 2 
g 



(IK 



(lUlj' 



(35) 



(36) 



In this case, the instability is dissipative since the growth 
rate is proportional to 77. 

If oj 2 + (i(k 2 + oj 2 ) > 0, then there should be no in- 
stability in the ideal magnetohydrodynamics (see condition 
(l27b ). Indeed, expressions d35l > and d36b yield Re7 = in 
the limit lo^ — > 0. However, in a dissipative MHD, the in- 
(30) stability can occur even if oj 2 + (i(k 2 + ui 2 ) > 0. In this 
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case, the first mode (non-oscillatory) is unstable alone if (we assume that the unperturbed disk is convectively stable 



uig + fin < 0, but two oscillatory modes are stable since us j 
should be positive. On the contrary, if uj 2 + [i(k 2 + wj) > 
but Wj < 0, then two oscillatory modes are unstable, but 
the first mode should be stable since ui 2 + /in 2 is positive. If 
lo 2 + h(k 2 + uj 2 ) < 0, then mode 3 is unstable with a very 

large growth rate 73 = \J\u 2 + h(k 2 + uj 2 )\ but mode 2 
is rapidly decaying with the decay rate approximately equal 
to 73. The first mode can also be unstable in this case if 
lo 2 + fin 2 > but its growth rate is small since it is propor- 
tional to the small dissipative frequency uj n . 

If u n is greater than other characteristic frequencies, 
then we have from Eqs. (f33t-(f34b 



7i 



72,3 



2uj n 



(37) 



Mode 1 is always stable but modes 2 and 3 can be unsta- 
ble. If lo 2 + /iK 2 > 0, then modes 2 and 3 are oscilla- 
tory. The instability of these modes occurs if uj j < 0. In 
the opposite case uj 2 > 0, oscillatory modes do not arise. 
If co 2 + fiK 2 < 0, then oscillatory modes become non- 
oscillatory, and one of these modes is unstable. The case 
of large ui v corresponds to small magnetic Reynolds num- 
ber and is of particular interest for protostellar disks. The 
instability of oscillatory modes can operate even in the dead 
zone of protoplanetary disks where the magnetic Reynolds 
number is small. 

3.2 Stability in a strong magnetic field with a e > 1 

Let us consider the stability of a strongly magnetized plasma 
with a e > 1. In this case, the growth rate is described by 
Eq. $1% with complex coefficients. The roots of Eq. $1% 
can be calculated by making use of general expressions 
(T3~3ll-(r3"4Tl for the roots of a cubic equation. However, these 
expressions are rather cumbersome and inconvenient for 
analysis. Therefore, we consider in detail the growth rate 
in the case when the frequency associated to electric cur- 
rents ojj is lower than the angular velocity £1 or character- 
istic buoyancy frequency ui g . The dissipative frequency uj v 
can be high and comparable to (or even higher than) other 
characteristic frequencies. This case is of particular interest 
for the dead zones of protostellar disks where the conduc- 
tivity is extremely low, and the magnetic Reynolds number 
can be relatively small. We can rewrite Eq. ( fT9l ) as 



.-' .2 ....2, 7AK; 



(7 +< + mO + 



= 0. 



(38) 



7 + IUJ 

This shape is more convenient to calculate the oscillatory 
roots (which can generally be unstable) by making use of 
a perturbation procedure. Since the last term on the l.h.s. is 
proportional to the square of a low frequency u>j, it can be 
considered as a small perturbation. Therefore, the solution 

of Eq. d38l can be represented as a power series of uj 2 : 7 

M 



where 7^°) does not depend on uj 2 and 7W 



-yW +7W 4- 

is linear in uj 2 . The equation of the zeroth order yields 



(39) 



and uj 2 + [ik 2 > 0). The roots are imaginary in the zeroth 
approximation, and there is no instability in the absence of 
electric currents. The correction of the first order is 
1 



7 



(i) _ 



(40) 



2 7 (°) + iuj ' 

Splitting this equation into real and imaginary parts, we 
have for the growth rate 



Re7 = — -- 



UJri/IUJj 



(Wl± 



+ [IK 2 ) 2 + UJ 2 



(41) 



Like the case of a weakly magnetized disk, the instability 
arises only if the magnetic field satisfies condition ((28). The 
growth rate depends on the wavelength of perturbations and 
can be essentially different for different k. If the wavelength 
A = 27r/fc is sufficiently short such as uj n > 57, then the 
growth rate is approximately given by 



Re7 



[1LU 2 



(42) 



One should note that the last expression (obtained for a case 
of strongly magnetized electrons) is similar to the one found 
in section 3.1 [see Eq.d37bl. It is assumed here that a dissi- 
pative frequency u> n is greater than all other characteristic 
frequencies. The order of magnitude estimate of Re7 is 



Re7 - 4.3 x 10 



-12 



-,1/2 



(43) 



where B, 



ip2 



nuT 2 

B^/IOOG and x e _ 12 = x e /lQ- 12 . The con- 



(44) 



12 



dition u>n > £1 is equivalent to 

A12 > 0.66T 2 1/2 P y 

where P yr is the rotation period in years and A12 = A/10 
cm. 

The imaginary part of 7 that determines the frequency 
of perturbations is approximately given by expression of 
the zeroth order d40| i. This expression describes buoyancy 
waves modified by differential rotation. It is well known that 
the buoyancy waves can be unstable and are responsible for 
convection if uj 2 < that is not likely to be satisfied in pro- 
tostellar disks. The instability considered in this section is 
the instability of buoyant waves as well and, in fact, is an 
oscillatory modification of convection. In stellar hydrody- 
namics, an oscillatory convection is often called semicon- 
vection and can be caused, for example, by a gradient of the 
chemical composition, As it is seen from our consideration, 
the distribution of the toroidal field can also be the reason 
of semiconvection in protostellar disks. 

In the opposite case, cu v << il, we have from Eq. d42l 
for the growth rate 

Re 7 «4 J ^- (45) 



2 co 2 



[IK 



The instability can occur for perturbations with u> n < £1 
as well, however, the growth rate is low in this case. The 
frequency of such weakly unstable perturbations is given 
by Eq. ©. 
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4 Summary and discussion 
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This paper examines the instability of weakly ionized pro- 
tostellar disks threaded by an external magnetic field. The 
conductivity of such disks is low, and dissipative effects can 
play an important role in the evolution of perturbations. This 
concerns particularly the midplane that is shielded from cos- 
mic rays to such extent that the MRI does not occur even 
under the most favorable conditions (Turner et al. 2006). 
Since the problem of stability is rather cumbersome with 
taking account of dissipative effects, we have considered a 
special case of perturbations with the wavevector perpendic- 
ular to the magnetic field. Such perturbations are not subject 
to the magnetorotational instability, because its growth rate 
is proportional to (fe • B) and is vanishing for the considered 
perturbations. 

It turns out that dissipative process can alter drastically 
the stability properties of protostellar disks. Apart from the 
instabilities that are typical for non-dissipative differentially 
rotating disks, new instabilities can occur that are deter- 
mined by dissipative processes. In the simplest case con- 
sidered in this paper, condition of the dissipative instability 
(l30l is entirely determined by the radial dependence of the 
azimuthal field and can be represented as 



din B u 



d In s 



> 1. 



(46) 



The instability arises if the magnetic field decreases with s 
faster than 1/s or increases more rapidly than s. Note that 
the condition of the considered instability does not depend 
on the magnetic diffusivity, and can be satisfied in the lim- 
iting cases of very high (rj — ► 0) and very low conductivity 
(?/ — > oo). However, the growth rate of instability depends 
sensitively on conductivity and is proportional to rj and 1 /rj 
in high- and low-conductivity limits, respectively. 

The condition of instability d37| i does not depend di- 
rectly on the rotation law. However, differential rotation can 
influence this condition indirectly because the radial profile 
of B v depends on fi(s). If stretching of the azimuthal field 
lines in the basic state is balanced by ohmic dissipation, then 
we approximately have from the induction equation 



r)AB v ~ sn'B s , 
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or 



B u 



s 3 n'B a 



(48) 



If rotation is Keplerian and Q oc s~ 3 / 2 , then the radial de- 
pendence of B v is given by B v oc s 1 / 2 B s /rj. Therefore, 
the considered instability may occur in the disk if the ra- 
tio B s /rj decreases with s faster than s -3 ' 2 . Likely, this 
condition can often be fulfilled in protostellar disks, par- 
ticularly if the radial field component decreases as a dipole 
(oc s~ 3 ). Therefore, the considered instability can generate 
turbulence in regions with a very low conductivity including 
the dead zones which likely exist in protostellar disks. 
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